Dynamics of Brans-Dicke cosmology with varying mass fermions 
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In this paper, the cosmological dynamics of Brans-Dicke theory in which there are fermions with 
a coupling to BD scalar field as well as a self-interaction potential is investigated. The conditions 
that there exists a solution which is stable and represents a late-time accelerated expansion of the 
universe are found. The variable mass of fermions can not vanish exactly during the evolution of the 
universe once it exists initially. It is shown that the late-time acceleration depends completely on 
the self-interaction of the fermion field if our investigation is restricted to the theory with positive 
BD parameter u). Provided a negative uj is allowed, there will be another two class of stable solutions 
describing late-time accelerated expansion of the universe. 

PACS numbers: 98.80.-k, 95.36.+X 



I. INTRODUCTION 

Strong evidences from the current cosmological ob- 
servations such as supernovae type Ia(SNeIa)[l], cos- 
mic microwave background (CMB)[23 and large scale 
structure(LSS)f3'] converge upon the fact that the uni- 
verse is spatially flat and there exists exotic component, 
dubbed dark energy, which drives the speed-up expan- 
sion of the universe. Many scenarios have been proposed 
to explain the acceleration in the framework of general 
relativity (GR). The preferred and simplest candidate for 
dark energy is the Einstein's cosmological constant which 
can fit the observations well. However, it suffers from 
the so-called fine-tuning problem and coincidence prob- 
lem. On the other hand, the observations are not yet 
able to confirm that dark energy is indeed a constant. 
Actually, many dynamical models of dark energy have 
been studied extensively during the past over tenyears, 
such as quintessence [J|, phantom [5], quintom [6| (see 
also Ref.[7i] for a detailed review), tachyon j8,], Chaplygin 
gas [3|3, etc (see Ref.[li| for a comprehensive review on 
dark energy models). 

GR is a well tested theory in solar system, but it is 
poorly tested on cosmic scales. Therefore, it is interest- 
ing to ask whether the cosmic acceleration is caused by 
a modified theory of gravity. As is pointed out, the dark 
energy problem may be essentially an issue of quantum 
gravity [I4I. Although a complete theory of quantum 
gravity has not been established, some valuable ideas 
are thought to be the features of the theory of quan- 
tum gravity, for example, extra dimensions. Many fun- 
damental theories would induce scalar fields in the usual 
4-dimensional space-time, and these scalar fields are in 
general non-minimal coupled. It is truly remarkable to 
find that a candidate scalar field of the desired nature 



is provided by some fundamental theor}0. On the other 
hand, in recent years, the interest in scalar-tensor theo- 
ries of gravity as viable alternatives to GR is renewed. In 
particular, some authors have resorted to scalar-tensor 
theory in order to explain the present accelerating ex- 
pansion of the universe [l^ [ij] . As the simplest example 
and prototype of scalar-tensor theory of gravity, Brans- 
Dicke theory (BD) formulates the gravitational phenom- 
ena through the interplay between the metric tensor and 
a scalar field </) that controls the intensity of the gravi- 
tational constant G, while the coupling between (j) and 
matter is absent. Although there is no direct support 
from fundamental theories, it is hard to deny that it ap- 
pears to provide a small window through which one can 
look into phenomenological aspects of fundamental the- 
ories. 

Recently, the possibility that fermion fields as gravita- 
tional sources could be responsible for accelerated peri- 
ods during the expansion of the universe is considered 
in the framework of GR[ia| as well as in BD theory 
of gravity[la|. In the context of cosmology, due to the 
isotropy and homogeneity of the geometry of the uni- 
verse, the fermion field configurations should be constant 
in space. It is shown that the fermion field could behave 
as dark energy for late-time universe. In this work, we 
would like to investigate the dynamics of BD cosmology 
in which there exists fermions with coupling to BD scalar 
field as well as a self-interaction potential, and the main 
focus is looking for the conditions under which there is a 
solution that is stable and represents a late-time acceler- 
ated expansion of the universe. 

This paper is organized as follows: In section |TT1 we 
briefly introduce the tetrad formalism of fermions with 
variable mass in Brans-Dicke theory of gravity. Then, 
the cosmological dynamics of the model is investigated in 
detail in section Hill Finally, in section ITVl we summarize 
our results and give some discussions. 
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^ It has been shown that only the generahzed version of the Chap- 
lygin gas model is compatible with the observations (see e.g. Ref. 



^ Note that string theory is a quite special type of scalar-tensor 
theory of gravity, and unless one introduces non-perturbative 
effects, the theory has a runaway problem. 



II. FERMIONS WITH VARIABLE MASS IN 
BRANS-DICKE THEORY 

We start from the Dirac Lagrangian density in 
Minkowski space-time 

Ld='i^ [^7"9aV - {dai>)l''^] - m^V - V, (1) 

where the index a = 0, 1, 2, 3, 7° is the Dirac-Pauh matri- 
ces, m is the fermion mass, f/i = ij}''^^^ denotes the adjoint 
spinor field and V , an exchisive function of ip and t/i, is 
the potential representing a fermion self-interaction. 

When gravity is taken into account, the generally co- 
variant Dirac Lagrangian becomes 

^D - ^ [^r'^i^^V - {D^i^)T^'^] - mV^V - V, (2) 



where R is the scalar curvature and is the Brans-Dicke 
scalar field. ifP' in the non-minimal coupling term acts 
as an effective gravitational constant Geff in such a way 
that G~^ — 27r0^/w. We assume that, just as the Brans- 
Dicke scalar field adjusts the strength of gravity G, 
the fermion mass m is replaced by /3(/), with /3 being a 
dimensionless coupling constant. Then the Lagrangian 
for the spinor field with variable mass can be written as 



ii3 = 2 [^r^i^^v- - (^p^)r''V] - mi^ - V, (6) 



So we will take the following total action 



S - / d^x,/^{LD + Lg + Lm), (7) 



where F^ — e^7° is generalized Dirac-Pauli matrices, I?^ 
denotes covariant derivatives, which is defined by 



D^^^d^ 



.9m^ 



[r;;, - eUd.ei)] 7^7' 



(3) 



where V^^ and e[^ denote the Christoffel symbol and the 
tetrad, respectively. The metric tensor g^^, satisfies the 
relation 



where Lm denotes the Lagrangian of matter fields. 

Varying action ([T]) with respect to the spinor field, we 
obtain the Dirac equations for the spinor field and its 
adjoint field 



ir'^D.i; - p<t>ij 



dV 



= 0, 



(8) 



9tJ.u 



^fj,^i/Vab: 



(4) 



where rjab is the Minkowski metric tensor. 

In the canonical frame, the Lagrangian for the Brans- 
Dicke theory can be written as 



8lu' 



LG = -i-R + -g^''y^c^y, 



(5) 



iD^ipV^" + P(j)i(j 



dV 

9V 



(9) 



Similarly, the Einstein equations for Brans-Dicke gravi- 
tational dynamics can also be derived from action ([7]), 



4w V^''" 2^^"^ 



2uj 



(g^,0n0 + 5^.(^0)' - 0V^v,0 - v^0V, 



rp4, 



T,f 



rpTi 



(10) 



r 



where D = g^'^V ^Vv is a covariant D'Lambertian for a 
scalar field and 



1 



T;t=V^0V,<^--5^,(V(/.) 



(11) 



- gfi^Lo- (12) 

The energy momentum tensor for the matter fields is 
defined as usual, T^ = -^ gl„ ■ In cosmology, it 
is often be expressed as the form of perfect fiuid 

TJZ^ iPrri+P7n)Uf_,U^ -Pmgfiiy, (13) 

where p„j and p„j are energy density and pressure of the 
matter, respectively. The four velocity vector m^ nor- 
malized as u^u'^ = 1. Finally, we vary action ([7]) with 



respect to 0, obtaining the equation of motion for the 
Brans-Dicke scalar field 



Dcl)- 



(t>R 

Aijj 



■ /JVV' = 0. 



(14) 



HI. 



DYNAMICS OF FERMION FIELD IN 
BRANS-DICKE COSMOLOGY 



For a spatially flat homogeneous and isotropic uni- 
verse, the space-time interval is written as usual as 



ds^ = dt^ - a^[t){dx^ + dy^ + dz^), 



(15) 



where a{t) is the cosmic scale factor. Therefore, the com- 
ponents of tetrad becomes 



-0 ~ "0 ' 






1 



and the Dirac matrices turn out to be 

a{t) 
from which the covariant derivatives is obtained 



I?o = ao, A = 5, - -a(i)7V, 



(16) 



(17) 



(18) 



where and thereafter overdot denotes derivative with re- 
spect to cosmic time t. 

Then, the Dirac equations ^ and © reacQ 

7^ + ^ffV + z/37>V + i7°^=0, (19) 



i^ + ^Hi,- i^fi,i, - ^7°!^ = 0, (20) 

and the equations of motion for the metric and the Brans- 
Dicke scalar field become 



4a; V a 



4a; zoj 2 



1 :^ 



(21) 



1 :^ 



^ 2- + i/^ +_iJ00+ 00+ 0^ = _p^_ 0^_p^ 



2a; 2a; 



(22) 



3if0 + /3^V-7^ f-+^')0-O (23) 
2a; V a 



where H = aja is Hubble parameter and 






(24) 



(25) 



Note that, among the equations (|2ip - (P5|) . only two of 
them are independent and the rest one can be derived 
from Bianchi identity. From Eqs. ([2T1) - (|23|) . we obtain 



a 
a 




4a; 

302 



+ 2{1+UJ)H^-:^{V+Pra)- 



and 



2ujp^ptp 

3(/)2 

(27) 



^ + 3if- + -^ = 



2a; 



02 (3 + 2a;)02 
dVtp 



-3 



^ 2'd^ 



(28) 



We shall consider a pressureless dust matter field, that 
is, Pm = 0. Then, from conservation equation for the 
energy density of the matter constituent, pm + •iH^pm + 
Pm) = 0, we can obtain that pm = PmoCL~^ ■ In order 
to analyze the cosmological dynamics, we should select a 
specific potential V . According to the Pauli-Fierz theo- 
rem, V can be an exclusive function of the scalar invariant 
(■0'(/')^[l5|. In this work, we choose 



V = X{il)il)f 



(29) 



as in Ref. [16| | , where A is a constant and a is a real num- 
ber. The bilinear ■0V' is scalar, hereafter we call it /. 
From Dirac equations ^9]) and (|20p . it is not difficult to 
find that / = foa^^, where /o is present value of /. 
Let us define some new dimensionless variables y — 

tions of X = In a. Therefore, Friedmann equation ([25]) 
becomes a constraint equation 



2a; 2 
1 = — y 
3 ^ 



2y + i}v + ^f + ^ 



(30) 



After a lengthy but straightforward calculations, from 
Eas. (P5)) - P5)) . it is shown that the deceleration parameter 



aH^ 



can be indicated as 



9=1 



2a) 



n 



-y 



/ 



3 + 2a; 



a;(4-6a) 

"m, „ "l^, (31) 



2a; 



where the evolution of y satisfies 

3 



y' = yq-'2y -"^y 



6-f 4a; 



[nra + {4.-6a)nv], (32) 



where the prime denotes derivative with respect to x. 

On the other hand, according to the definitions of ftf 
and rim, as well as the evolutions of / and pm, it is easy 
to find that 



n'=-il + y-2q)nf, 



(33) 



^2 _ 2a; 02 



H'' = ^^-2H^ + ^{V + /3(f>^^, + p,n), (26) 



nL 



{l + 2y- 2q) n. 



(34) 



^ Just as mentioned in the introduction, for the isotropic and ho- 
mogeneous universe, the fermion field is an exclusive function of 
time. 



The equations ([32)) - ([34)) constitute an autonomous sys- 
tem and we shall use phase space method to investigate 
the properties of its solutions qualitatively. The critical 
points of the system can be obtained by setting il'r — 0, 
il'^ = and y' = 0. In Table HI all the critical points 



of the system are listed. The coordinate of the point D 
is independent of w and a, while the points Ji± and C 
is independent of a. As for the coordinates of the other 
two points A and E, they rely on both of the two pa- 
rameters w and a. Noting that merging occurs for the 
points D and E for a — ^. All the critical points except 
C are associated with the dust matter vanishing universe. 
From Eas. ([55)) and (|M)) . it is obvious that plane y — ftm 
(flf = 0) and plane y — flf (flm = 0) are two invari- 
ant sub-manifold of the 3-dimensional phase space. This 
means that there is no orbit in the phase space for which 
rim or ri/ can be exactly zero. That is to say, if the 
initial value of / = ipi^ or coupling constant /? in the 
definition of flf is set to be zero, i.e., the fermion field 
is initially massless, it will remain massless. The values 
of fly and deceleration parameter at the different critical 
points flvc and Qc are also listed in Table HI 

What is more interesting is the stable solution which 
represents that the universe will undergoes a late-time 
accelerated expansion phase. To investigate the stability 
of these critical points, we can write the variables near 
these points (fife, ^mc, yc) in the form y = yc + t]!, il/ = 
^fc + m and ^rn = ^mc + V3 with rj, {i = 1,2,3), the 
perturbations of the variables around the critical points 
to the first order. This leads to the following equation 



U' = A • U, 
where the 3-colunm vector U — {rjiY 



(35) 



sents the perturbations of the variables and A is a con- 
stant 3x3 matrix. For stability we require all 3 eigen- 
values of A to be negative. 

The eigenvalues of the linear perturbation matrix A for 
the critical points and the stability properties are sum- 
marized in table [ITl while the conditions for a stable and 
late-time acceleration solution are listed in table IIIIl It 
is found that all the critical points but B_ can be stable 
for some values of parameters lo and a. However, the 
stable solutions presented by D and E can not describe 
a universe that undergoes a late-time acceleration. Fur- 
thermore, it is worth noting that only point A allows BD 
parameter w > for a stable and late-time accelerated 
solution. 

As is well known, in the limit case cj — > oo, the 
standard Einstein's GR is recovered from the BD the- 
ory, except for some special cases, for example, in which 
the trace of the matter energy-momentum tensor van- 
ishes [i7|. And the solar-system experiments constraint 
severely that the value of u) obeys cj > 4 x 10'' jl8| . 
With this understanding, the late-time acceleration of ex- 
pansion of the universe will be entirely dependent upon 
the self-interaction of the fermion field as is illustrated 
in figure [1] where the BD parameter is chosen to be 
Lu — b X 10^, and initial values of y, 0/ and Vim is 
set at present (a = 1) to be -0.00001, 0.05 and 0.3, 
respectively. We recall that in the standard FRW uni- 
verse containing only single barotropic fiuid in the frame- 
work of GR, the deceleration parameter q and the equa- 
tion of state of barotropic fluid w = p/ p have a relation 




FIG. 1. The evolution of deceleration parameter q with re- 
spect to In a for different values of a, where the BD parameter 
is chosen to be a; = 5 x 10^, and initial values of y, Qf and 
Q,m is set at present (a — 1) to be —0.00001, 0.05 and 0.3, 
respectively. 



2q = 1 + 3w. Therefore, from figure [1] we can find that 
the equation of state of effective dark energy is able to 
cross the phantom divide w = — 1 for some values of a 
and the big rip may also occur in the future. In figure [21 
we plot the evolution trajectories of relative energy densi- 
12 3 repre- ties 51^, Qf and Qy- It is not surprising that the energy 



density of the fermion field will eventually overcome the 
energy density of the matter field and this occurs roughly 
at the time the universe going into an accelerated expan- 
sion phase, while the energy density contributed by the 
variable mass of fermion field keeps similar pace with that 
of the matter field. 

On the other hand, if we allow the BD parameter uj 
to be less than zero on the cosmological scales, there 
will exist another two class of stable late-time acceler- 
ated expansion solutions which are represented by criti- 
cal points B+ and C and illustrated in figure [3] and figure 
m respectively. These solutions do not determined by 
the self-interaction of the fermion field, because in these 
two class of solutions J7y converges towards zero and the 
parameter a > 1/2. The solutions illustrated in figure 
13] denote a phantom field dominated universe which will 
suffer from the fate of big rip. However, the solutions 
illustrated in figure |4] will avoid this violent doom, al- 
though the equation of state of effective dark energy un- 
der certain conditions may become less than —1 within 
some period of time. More interestingly, it is shown in 
figure [5] that, in this class of solutions, the relative energy 
density of dust matter will not vanish but converge to a 
reasonable finite value. 



IV. CONCLUSIONS 

Brans-Dicke theory is a natural alternative and a sim- 
ple extension of general relativity. In this work, we have 



Points 


Coordinate {Qfc,0,mc,yc) 


nvc 


9c 


A 


(0,0,,,:^^2a_) 


{3+2bj)l^5+6uj+6a'-'uj-6a{l+2uj)') 
3(l-2(a-l)t^)^ 


2(ci-l)cj-l 




(o,o,|,(i±Vi + 2f)) 







B± 


3 + 2aj±V9+6i<j 


C 


V*^' 6(l+uy^ '2^+2J 





2+uj 
2+2ij 


D 


(1,0,0) 





1 
2 


E 


(6a - 2 - 12uj + 36aui - 24a'^uj, 0, 3 - 6a) 


3{l-2a)(3 + 2uj) 


2-3q 



TABLE I. The critical points of the autonomous system Eas. p2|l - (|34|l . The values of ^v and deceleration parameter at the 
different critical points flvc and qc are also listed. 



Points 


Eigenvalues 


Stability 


A 


5 + 6tJ + 6a^ai - 6a(l + 2aj) 

-l + 2{-l + a)ui 

1 + 6cj + 12a^Lj - 3q(1 + 6ti)) 

-1 + 2(-l + Q!)W 


1 / -5 + 6a -1 + 3q \ 
attractor for a < — & uj < or l<j > ; 

2 V 6 - 12a + 6a2 6 ~ ISa + 12q2 y 

/2 _1 + 3q _5 + 6a \ 

rcpcllcr for - < a < 1 & < l^ < 

V 3 6 - 18a + 12a2 6 - 12a + Ga^ J 

f -l+3a -6+ea \ 

or a > 1 & U. < or u> > ; 

V 6 - 18a + 12a2 6-12a+6a2y 

saddle for other cases. 








B+ 


3 + Sw ± VS + 6w 


/I -5 + 6a ^/l 4 \ 

attraetor for a < - & <uj<Oora>-& < l^ < ] ; 

V 2 6 - 12a + 6a2 / V 2 3 / 


3 (3 + 2w ± V9 + 6a;j 

2aj 


/ -5 + 6q \ 

repcller for (a < 1 & c^ > 0) or a > 1 & < l<j < ; 

V 6-12a+6a2; 




3 - 6(a - l)aj ± V9 + 6aj 








B_ 


3 3 + 6uj - V9 + Giu 

rcpcllcr for oj > & a < ; 

2 6u; 

saddle for other cases. 


C 


9 fin 4+3t>J 1 


1 4 


" ua, 2+2i^ ' 2+2uj 


2 3 
saddle for other cases. 


D 


-|,0,3-6q 


attractor for a > ^; saddle for a < ^. 


E 


6a - 3, 3q - 2 ± ^ 

whcro A = 8 + 24cj - 8o,(4 + 15u) 


1 S - 32a + 33a2 3a - 1 
attractor for a < — & < uj < ; 

2 96a3 - 192a2 + i20a - 24 6 - 18a + 12a2 

2 8 - 32a + 33a2 3a - 1 
rcpcllcr for - < a < 1 <^ < i^ < 

3 96a3 - 192q2 + 120a - 24 6 - 18a + 12a2 

3a - 1 8 - 32a + 33a2 

or a > 1 &: < uj < ; 

6 - 18a + 12a2 9Ga^ - 192a2 + 120q - 24 

saddle for other cases. 



TABLE II. The eigenvalues of the linear perturbation matrix and stabihty for aU the critical points. 



Points 



Conditions for a stable and late-time acceleration solution 



A 



a < - &c {u < -\uj > - 

3 6 - 12a + 6a2 2 - 8a + 6a^ 

1 „ 9 

or a — — &u;< 

3 8 





or 


1 
3 


< 


1 „ 1 -5 + 6q 

a < - & < ui < - 

2 2-8q + 6q2 6-12Q + 6a2 


B+ 


-f <ui <Oba> ^i; 2_ 


C 


-2<u;<-|&a>i 


D 


False 


E 


False 



TABLE III. The conditions for the critical points representing a stable and late-time accelerated solution. 





FIG. 2. The evolution of relative energy densities with respect 
to In a, where a and uj is respectively chosen to be —0.05 and 
5 X lO** and the initial (present) values of j/, 57/ and Q,m is set 
by -0.00001, 0.05 and 0.3. 



FIG. 4. The evolution trajectories of deceleration parameter 
q with respect to In a for different initial (present) values of 
y, where a and uj is respectively chosen to be 0.9 and —1.95. 




FIG. 3. The evolution trajectories of deceleration parameter 
q with respect to In a for different initial (present) values of 
y, where a and uj is respectively chosen to be 1 and —1.25. 
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FIG. 5. The evolution trajectories of relative energy density of 
dust matter with respect to In a for different initial (present) 
values of y, where q and lo is respectively chosen to be 0.9 
and -1.95. 



studied the cosmological dynamics of Brans-Dicke the- 
ory in which there are fermions with a coupling to BD 
scalar field as well as a self-interaction potential. Because 
of the isotropy and homogeneity of the geometry of the 
universe, the fermion field is independent of space and 
becomes an exclusive function of time. The conditions 
that there exist solutions which are stable and describe a 
late-time accelerated expansion of the universe are found. 
The variable mass of fermions can not vanish exactly dur- 
ing the evolution of the universe once it obtains initially. 
It is shown that the late-time acceleration depends com- 
pletely on the self-interaction of the fermion field if our 
investigation is restricted to the theory with positive BD 
parameter uj. Provided a negative a; is allowed, there 
will be another two class of stable solutions describing 
late-time accelerated expansion of the universe which are 
independent of the self-interaction of fermion field. 



The strongest constraint to date on the Brans-Dicke 
theory of gravity has been put on the solar system 
scale [1^, however, it should be pointed out that such con- 
straint does not necessarily apply on scales much larger 
than those of the measurements, and on epochs much 
different from the present. The solar system experiments 
only probe scales in gravitational equilibrium, where the 
background expansion of the universe is negligible, there- 
fore, they would not reveal spatial or time variation of the 
gravitational constant on larger scales. Through a certain 
mechanism, the BD parameter w becomes variable from 
small scales to large scales. It is even not excluded that 
on the largest scales the gravitational constant becomes 
effectively negative. Some indications for a negative w 
(or effectively negative gravitational constant) have been 
obtained [1^. Therefore, it is interesting to investigate 
such a model presented here in the context of a detailed 



viable gravitational theory with variable effective w, for 
example, the abnormally weighting energy theory [iQ"!. 
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